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A slip velocity method for two-dimensional incompressible turbulent 
boundary layers was presented in reference 1. The inner part of the boundary 
layer was characterized by a law of the wall and law of the wake, and the 
outer part was characterized by an arbitrary eddy viscosity model. In the 
present study for compressible flews, only a law of the wall is considered. 
The problem of two-dimensional compressible flow is treated first? then the 
extension to three-dimensional flow is addressed. 

Two-Dimens ional Compressible Flow 

Basic Equations 

The governing eguations for compressible boundary layer flow are 
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where x and y are the normal and tangential coordinates, u and v are 
the respective velocity components, p is the density, p Q and U are the 
density and fluid speed at the outer edge of the boundary layer, and p is the 
dynamic viscosity. The sum of the dynamic viscosity and the turbulent eddy 
viscosity p^. is defined as 

p + p fc = K(x,y ) p e U6* (3) 

★ 

where K is a general nondimens ional function of x and y, and 6^ is the 
incompressible displacement thickness. 

In this treatment, the defect stream function of Clauser (ref. 2) is 
used. This function is defined as 
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The prime denotes differentiation with respect to r\, and the shear stress 
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velocity u is defined as 
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where x and p are the wall shear stress and density. The boundary layer 
w w 

thickness parameter A is defined as 
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Partial derivatives with respect to x and y are of the form 
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where the dot denotes differentiation with respect to £. 

Law of the Wall 

It is assumed that the flow is adiabatic, and the gas is calorically 
perfect. The law of the wall for this flow is obtained with a treatment 
similar to that of Van Driest (ref. 3) as 
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where y is the ratio of specific heats, tc and B are the von Karman 
constants, and a QW and p gw are the adiabatic-wall speed of sound and dynamic 
viscosity, respectively. The quantity r is the recovery factor, which is 
typically evaluated as , where Pr is the Prandtl number. The equation 

for the density for this flow is 
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where M is the edqe Mach number. Since the ratio u /a . is small, the 
“ aw 

equation for u can be written as 
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Note that the lowest order term in this equation is the same as the law of the 
wall for incompressible flow. 

Hie nondimens ional shear stress velocity ratio is 


Y = U 


and where m(5) is defined as 
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The Clauser pressure gradient parameter p is defined as 
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where p is the pressure. The dot can be used to represent the qradient of 
x as well as £ since U is a function only of x = 5. 

An important relationship for the lowest order form of the stream 

a 

function f(£»n) and its derivative f is obtained from the law of the 
wall. Equation (8) can be written in terms of n and K as 
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where the Reynolds number based on the incompressible displacement thickness 
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6., the wall properties p and p , and the edge velocity U is 
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The stream function f can be expanded in terms of the small parameter y as 

f = f + iff ,+ . . . 
o 1 


With this expansion and equation (4), equation (12) can be written to lowest 
order as 
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This equation pertains throughout the inner region of the boundary layer. 
Governing Equation for f 0 

The governing equation for f Q is obtained from equation (2), the 
tangential momentum equation. To establish the equation for f Q , the mass flux 
components pu and pv and the partial derivatives of u with respect to x and y 
must be expressed in terms of f Q . With equations (4) and (7), the u component 
of velocity and the density can be expressed to first order in y as 
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With these equations the mass flux pu can be approximated as 


pu = p u {1 + y ( 2 ^ -1) f '} 
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From equations (6) and (14) the derivative 3u/3x is obtained as 
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where the nondimensional tangential coordinate s is defined as 
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The zero-order approximation for the mass flux component pv is obtained from 
equations (1), (6), and (15) as 
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Finally, the normal derivative of u is 
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With equations (15), (16), (17), and (18), the tangential momentum equation 
can be written as 
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It should be noted that this stream function treatment is patterned after 
that of Mellor and Gibson (ref. 4) for incompressible equilibrium boundary 
layers. In particular, it is patterned after the lowest-order treatment of 
Mellor and Gibson. 

The three boundary conditions for f Q involve the values of f Q at the wall 
and in the free stream and the value of the shear stress at the wall. Since 
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With boundary condition (22), this equation can he integrated across the 
boundary layer to obtain 
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With this value, the first integral of equation (19) for arbitrary ri is 
written as 
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Note that no assunptions have been made which limit the arbitrariness of the 
nondimens ional viscosity coefficient K. 

Match Point Location 

The match point divides the outer and the inner reqions of the boundary 

layer. In the outer part of the boundary layer the viscosity coefficient K is 

arbitrary. In the inner region, the flow is governed by an empirical law 

of the wall. There is one point, the match point, at which both the arbitrary 

coefficient of viscosity and the law of the wall pertain. At this point the 

stream function f Q and its first three derivatives with respect to n are 

continuous. Note that in the parlance of asymptotic expansions, the "match 

point" would be properly termed a "patch point." 

It is assumed that the flew in the inner region and hence at the match 

point is essentially in equilibrium. The quantities f" and f 

o o 

in the inner region are evaluated with the law of the wall as 
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the mass flux component pv must vanish at the wall, the defect stream 
function f Q must also vanish at the wall: 



The free stream boundary condition is obtained from the definition of the 
incompressible displacement thickness as 
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With equations (4), (5), (9), and (18), the shear stress and shear stress at 
the wall can be written as 
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From these equations, it is seen that the shear stress boundary condition is 
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The qoverninq equation for f Q has a first inteqral. Equation (19) can be 


written as 
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Thus at the match point the governing eguation (23) can be written as 


Pp» v ,1 * 

— - — • 1 - (i + 2 — $£ 
p w ' <V °' m 



Let 


A= 



i 


Bf 


o,m 


and note that 


r m 


n m = A 
m A 


The governing equation is 
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Shear Stress 


The governing equation in the inner region can be written as 


If u and 3u/3y are known at sane point y, the shear stress can be determined 
as 
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Edge Velocity Determination 

If the total pressure at the edge of the boundary layer can be defined, 
it can be used to determine the edge velocity U. The tangential momentum 
equation at edge of the boundary layer is 
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where p fc is the total pressure at the boundary layer edge. 

A second approach is to define the boundary layer edge in tern© of the 
deviation of the total enthalpy from the freestream value. 

Solution Algorithm 

This solution process is for either an iterative solution of the 
viscous-inviscid problem or a marching solution of the viscous problem with 
the inviscid solution known. In either case, approximate values for the 
solution u(x,y), v(x,y), p(x,y), and p(x,y) are known. Also, the parameter u* 
is known. The six steps for one iteration are; 

(1) From the turbulence model compute the edge velocity U and the product 
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U6. 

(2) Calculate 
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(3) Compute iteratively 



11 


(4) Calculate 
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(5) Compute the new solution u(x,y), v(x,y), p(x,y), and p(x,y) 

usinq u and v as boundary conditions, 
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(6) Calculate 
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Note that the boundary conditions u and v can be applied at the model 

mm 


surface as a first approximation. 

The main difficulty with the present method is the need to specify the 
edge velocity U. For hypersonic flow and slender body flow U can be 
replaced with u^ . 


Cross-Flow Effects 

Initially, at least, it will be assumed that a small crossflow approach 
can be used. Cooke (ref. 5) has shown that the equation for the streamwise 
velocity component, measured relative to the velocity vector at the outer edqe 
of the boundary layer, for incompressible flow is independent of cross-flow 
effects and hence is similar to the equation for two-dimensional flow. It 
will be assumed that the same is true for compressible flow. 
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In this treatment , it is assumed that the turning angle near the wall is 
essentially linear with respect to distance from the wall so that it can be 
evaluated by extrapolation. The linear behavior of the turning angle near the 
wall is supported by the data of Johnston (ref. 6) and Van Den Berg, Elsenaar, 
Lindhout, and Wesseling (ref. 7). 

An additional arguement for linear extrapolation comes from the results 
of Mellor (ref. 8) for two-dimensional, incompressible, high Reynolds number 
flow and the three-dimensional extension of Goldberg and Reshotko(ref . 9). 

Both of these treatments show that, to lowest order, the inner layer flew is 
determined by viscous forces; pressure-gradient forces do not appear* Since 
turning is the result of the interplay of pressure and viscous forces within 
the boundary layer, it is reasonable to assume that turning is complete by the 
time the inner region is reached and only viscous forces remain. 

For purposes of illustration, let x,y, and z be Cartesian boundary layer 
coordinates with y the normal coordinate. Ihe respective velocity components 
are u, v, and w. The turning angle $ and tangential velocity u fcan 
are defined as 

, . -1 ,w. .2.2. 1/2 

<|> = tan (— ) , u fcan = (u + w ) 

Let y ! and y 2 be the two locations immediately above the match point y^. 
The match point turning angle is 


Vo ~ V, 


♦ - 


m , ‘ 1 


yi - y, 


m 


m y 2 - y x y 2 - V] 


and the law of the wall is expressed in terms of u fcan as 
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The velocity ccnponents and w^, which are needed as slip velocity boundary 

conditions, are obtained from $ and u. with the equations 
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Inplementation 

Consider an x,y, z Cartesian coordinate system fixed to the three- 
dimensional configuration with the x coordinate in the axial direction. Let 
the unit normal to the surface be 

n s 1 n x + 1 °y k n z 

The flow anqle $ will be measured in the tanqent plane (the plane normal 
to n ). The anqle 4 > will be measured from the line where a reference 
plane intersects the tanqent plane. The reference plane will be either the 
x-y plane or the x-z plane, depending upon whether the projection of the unit 
vector n on the y-z plane is more closely aligned with the unit 
vectors 3 or £ . 

The unit vector n in the direction of the projection of n on the y-z 
plane is 
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The conditions for choosing the reference planes are: 
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x-y Plane as Reference 

A vector in the x-y plane normal to n is 

e xy = n x k = l x R n x + 3 x k = in y - jn x 

The unit vector in the x-y plane normal to n is 
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The velocity component in the x-y plane tangent to the surface is 
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where the velocity vector is 
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The velocity component in the x direction tangent to the surface is 
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The crossflow angle 4> is defined as the angle in the tangent plane between 
the velocity vector and the intersection of the tangent plane with the x-y 
plane: 

tan -1 | - “V* - vn v n z *_ w(n x * n y> 
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Let n he the coordinate in the normal direction. The angle $ at 
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the surface is obtained by linear extrapolation of the values ^ at n^ and <j> 2 
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Now assume that the new "slip velocity" , the angle 6 , and the 
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"inner layer transpiration velocity" have been determined. The values u 

v, and w are needed as boundary conditions: 
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U = a cos $ , w = U. sin * 
c tan m c tan T m 


The velocity vector can be written as 
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The Cartesian velocity components are 
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x-z Plane as Reference 

Now use the x-z plane for reference. A vector in the x-z plane normal 
to n is 
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The velocity component in the x-z plane tangent to the surface is 


U = $ • 1 


un 


wn 


xz , 2 . 2.1/2 

(n + n ) ' 
x z 


, 2 2 . 1/2 

(n + n ) ' 
x z 


(un z - wn x ) 

T ~ 2 ~ 072 

(n + n ) ' 
x z 


The velocity conponent in the 1- direction tanqent to the surface is 


w = 0 . 1 - = 

s c 


Un xV V(n x + n z } + Wn z n y 


. 2 . 2 . 1/2 
(n + n_) 


The crossflow angle <|> is defined as the angle in the tangent plane between 
the velocity vector and the intersection of the tanqent plane with the x-z 
plane: 

♦ = tan (jj-) = tan { y ~ H 

S Z X 

The value of $ at the surface is determined by extrapolation as before. 

Now assume that d> and U. are known. It follows that 
Y m tan 


U = IT cos d> f W = U sin <b 
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Also the inner layer transpiration velocity V n is known. The velocity vector 
can be written as: 


=1 U+nV+i-W 
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The Cartesian velocity components are: 
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